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^ ■ Context. In the framework of the study of the Galactic metallicity gradient and its time evolution, we present new high-resolution 

' spectroscopic observations obtained with FLAMES and the fiber link to UVES at VLT of three open clusters (OCs) located within 

I ~7 kpc from the Galactic Center (GC): NGC 6192, NGC 6404, NGC 6583. We also present new orbit determination for all OCs with 

, Galactocentric distances (Rqc) ^8 kpc and metallicity from high-resolution spectroscopy. 

Aims. We aim to investigate the slope of the inner disk metallicity gradient as traced by OCs and at discussing its implication on the 
chemical evolution of our Galaxy. 

Methods. We have derived memberships of a group of evolved stars for each clusters, obtaining a sample of 4, 4, and 2 member stars 
in NGC 6192, NGC 6404, and NGC 6583, respectively. Using standard LTE analysis we derived stellar parameters and abundance 
ratios for the iron-peak elements Fe, Ni, Cr, and for the o'-elements Al, Mg, Si, Ti, Ca. We calculated the orbits of the OCs currently 
located within 8 kpc from the GC, and discuss their implication on the present-time radial location. 

Results. The average metallicities of the three clusters are all oversolar: [Fe/H]= +0.12 ± 0.04 (NGC 6192), +0.11 + 0.04 (NGC 
6404), +0.37 ± 0.03 (NGC 6583). They are in qualitative agreement with their Galactocentric distances, being all internal OCs, 
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QQ . and thus expected to be metal richer than the solar neighborhood. The abundance ratios of the other elements over iron [X/Fe] are 

' consistent with solar values. 

' Conclusions. The clusters we have analysed, together with other OC and Cepheid data, confirm a steep gradient in the inner disk, a 

signature of an evolutionary rate different than in the outer disk. 
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1. Introduction many strong observational constraints, as those provided by the 

■ ■ ' radial metallicity gradient and its evolution with time. 
Theoretical models of chemical evolution allow us to understand 

the processes involved in the formation and evolution of galax- In spite of many observational and theoretical efforts, vari- 

ies, and in particular of the Milky Way (MW). However they de- ous questions remain open even in our Galaxy, where plenty of 

pend on many valuables, which include: the star formation rate observations based on 8m-class telescopes of stellar populations 

(SFR), the initial mass function (IMF), radial inflows and out- of different ages ai-e available. We are fai- from obtaining a firm 

flows of primordial and/or pre-enriched gas, the star formation result about the shape of metallicity gradient in particular in the 

law and possible thresholds in the gas density for the formation inner disk and at large Galactocentric distances, and about its 

of stars. To decrease the number of degrees of freedom they need time evolution. 



* Based on observations obtained at the ESO Based on spectro- Contrasting results about the time evolution appear from 

scopic optical observations made with ESO telescope VLT, program: the comparison of abun dances of planetary nebulae of differ - 

ID: 083.D-0682, Title: The Galactic radial metallicity gradient in the ent ages (cf. Maciel et al. |2007i Stang hellini & Haywood [2010] l, 

inner disk Cepheids (e.g., Pedicelli et al. 120091 Andriewsky et al. I2002I I. 
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HII regions (e.g., Rood et al. 120071) and blue supergiant stars 
(e.g., Luck et al. [19891 Daflon & Cunha'2007). 

In this framework, the determination of the chemical compo- 
sition of stars belonging to open clusters (OCs) represents prob- 
ably the best tool to study the behaviour of iron and other ele- 
ments across the Galactic disk. OCs are indeed populous, coeval 
groups of stars with similar chemical composition; they are lo- 
cated all over the Galactic disk and span large ranges of ages and 
metallicities; furthermore, the estimate of their age and distance 
is affected by much smaller uncertainties than those of field stars. 

The first determination of the Galactic metallicity gradi- 
ent with the use of OC metallicities was performed by Janes 
( 11979 ). He based his study on a sample of OCs covering the 
Galactocentric distance (Rqc) range 8-14 kpc and on metallici- 
ties ([Fe/H]) estimated from photometry, deriving a gradient of 
-0.075+ 0.034 dex kpc Several subsequent studies were per- 
formed using a variety of techniques and larger samples of clus- 
ters ( Friel & Janes . 1 993 i Fri el [T995l Tw arog e t al. il997t Carraro 
et al. [19981 Friel et al. [20021 Chen et al. [2053l l. all agreeing on a 
negative slope between -0.07 and -0.1 dex kpc 

More recently, with the advent of high-resolution spectro- 
graphs on 8m-class telescopes, it has become possible to derive 
more secure [Fe/H] estimates, and also to extend metallicity de- 
terminations to OCs in the outer disk, beyond 15 kpc from the 
Galactic center (Carraro et al. 2004'; Carretta et al. 2004; Yong et 
al.|2005; Sestito et al. 2006 2008 ). These studies have confirmed 
the steep slope in the gradient for 7 kpc< Rqc <1 1 kpc and have 
shown that the distribution becomes flatter for Galactocentric 
distances (Rcc) above 11-12 kpc (cf. Magrini et al. 2009, here- 
after M09, but see Pancino et al. 2010 for a different view). 

However, in spite of the significant progresses so far 
achieved, a very interesting radial region, Rcc < 7 kpc, has been 
neglected during the past spectroscopic studies of OCs. This 
radial region is particularly interesting because from Cepheid 
metallicities, several authors (e.g., Andrievsky et al. 2002; Luck 
et al. 12006 ; Pedicelli et al. l2009l l. have shown that there is a fur- 
ther change of slope of the gradient. From ~4 to ~7-8 kpc (val- 
ues vary slightly in different papers) the gradient is much steeper 
than in the outer regions. The value of the slope of the very inner 
(^GC <7 kpc) gradient for iron varies from -0.13 to -0.15 dex 
kpc The questions are if this holds true also for OCs, if this 
change of slope is time-dependent (and OCs can tell us this be- 
cause they span a large age range), and where the metallicity 
stops increasing due to the presence of the Galactic bulge. So- 
far only one cluster at a distance from the GC smaller than 7 kpc 
has an available metallicity determination (Mil, Rgc=6.86 kpc, 
age=200 Myr, [Fe/H]= -h0.17 - Gonzales & Wallerstein 2000). 

Several studies and chemical evolution models have tried to 
explain the changes of slope in the abundance gradient, among 
them, e.g. i) the study of dynamical effect of corotation res- 
onance (located close to the solar Galactocentric distance) by 
Lepine et al. (2003) which is assumed to be the main cause of 
the formation of the bimodal radial distribution of metallicity ; 
ii) the model by M09 in which the inner disk is formed inside- 
out by the rapid dynamical collapse of the halo and the outer disk 
is formed at a lower rate, from a 'reservoir' of gas; Hi) the model 
by Colavitti et al. (I2009I I where the inside-out disk formation 
and the density threshold for star formation are both necessary 
ingredients to reproduce the change of slope; iv) the model by 
Fu et al. (120091 1 where the star formation efficiency is inversely 
proportional to the distance from the galaxy center 

The prediction of these models are also different for what 
concerns the inner gradient. They all reproduce one to a few 
changes of slope in radial gradient but at different Rcc ■ for exam- 



ple, the model adopted in the study by Lepine et al. (120031 ) has a 
flattening of the gradient for Rcc > Rq, while the model by M09 
produces a gradient with an outer plateau (Rcc > 11 - 12kpc), 
and a very steep gradient for Rcc < 7kpc. Finally, the best 
model by Colavitti et al. ( 120091 ) has a more pronounced gradi- 
ent in the inner disk for Rcc < 7 - 8kpc, while the models by 
Fu et al. (120091) show a change of slope in the inner Galaxy for 
Rcc<2- 3kpc. 

The aim of this paper is to extend the sample of OCs at 
Rcc <7 kpc, thus providing constraints suitable to disentangle 
the roles of different model parameters in inducing changes into 
the chemical gradient and its slope as, e.g., the importance the 
dynamical collapse of the halo in the inner regions with respect 
to the accretion due to infall from the intergalactic medium and 
from mergers. 

The present paper is structured as follows: in Sect.|2]the three 
OCs under analysis are presented. In Sect. [3] we describe the ob- 
servations and data reduction. In Sect. |4] we show the abundance 
analysis, while in Sect. 5 we present our results and compare 
them with the inner disk giant stars, whereas in in Sect. 6 we 
discuss their implication in the shape of the metallicity gradient 
and in the chemical evolution of the Galaxy. In Sect.|2]we present 
orbit calculations and their consequences for the metallicity gra- 
dient. Finally, in Sect. [8] we give our summary and conclusions. 



2. Target clusters 

Our sample includes three of the most internal open clusters 
(age < 1 Gyr) whose evolved stars can be studied with a ground- 
based 8m-class telescope. They are NGC 6192, NGC 6404, and 
NGC 6583. 

Several photometric studies of NGC 6192 are present in the 
literature. These studies derived, however, different reddening 
values and, consequently, discrepancies in ages and distances ex- 
ist. The reddening values in the literature range from E(B-V) = 
0.26 (Kilambi & Fitzgerald [19831 King [T987b to 0.68 (Kjeldsen 
& Frandsen 11991), with the most recent estimate by Paunzen et 
al. (200j|) giving E(B-V) = 0.54. The related estimates of age 
range from ~1 Gyr to 0.09 Gyr, and those for the distance to 
the Sun from ~1 kpc to 1.7 kpc. Paunzen et al. ( ,2003 ) estimated 
also the metallicity of the cluster: [Fe/H] = -0.10 + 0.09. Loktin 
et al. (120011 ) revised the original data of Kjeldsen & Frandsen 
(1991]), and determined a reddening of E(B-V) = 0.64, a dis- 
tance of 1.5 kpc and an age of 0.13 Gyr. Five stars of NGC 6192 
(Nos. 9, 45, 91, 96 and 137) have CORAVEL radial velocities 
in the narrow range -8.8 km s ' < Vr < -6.4 km s ', appearing 
thus to be red giant members of NGC 6192 (Claria et al. 120061 ). 
Claria et al. (2006) also estimated the age (0.18 Gyr), the dis- 
tance (1.5 kpc), and the DDO abundance index finding a cluster 
metallicity [Fe/H] = +0.29+0.06. 

NGC 6404 and NGC 6583 have been studied by Carraro et 
al. (2005) with CCD photometry. From isochrone fitting they 
derived the main parameters of the two clusters. They found that 
NGC 6404 is 0.5 Gyr old and located at a distance of 1.7 kpc 
from the Sun towards the Galactic Centre direction, while NGC 
6583 is 1.0 Gyr old and at 2.1 kpc. E(B-V) values are 0.92+0.05 
and 0.51+0.05 for NGC 6404 and NGC 6583, respectively. For 
both clusters, solar metallicity isochrones provide a reasonable 
fit across the whole CMDs. No radial velocities are available for 
the stars of these clusters. 

In Table 1 we report the cluster parameters and references. 
We also calculated the Galactocentric distance, Rgc, adopting a 
distance of the Sun from the Galactic centre of 8.5 kpc. 
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Fig. 1. A region of the normalized spectra of member stars in NGC 6192, NGC 6404, NGC 6583. 
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Table 1. Cluster parameters from the literature. 



Cluster 


Age 


[Fe/H] 


Rgc 


d 


E(B - V) 


Ref. 




(Gyr) 




(kpc) 


(kpc) 






NGC 6192 


0.18 


+0.29 


7.1 


1.5 


0.64 


a 


NGC 6404 


0.5 


solar 


6.8 


1.7 


0.92 


b 


NGC 6583 


1.0 


solar 


6.4 


2.1 


0.51 


b 



a) Claria et al. ( l2006t 

b) Carraro et al. ( 12005) 



3. Observations and data reduction 

The three OCs were all observed with the multi-object instru- 
ment FLAMES on VLT/UT2 (ESQ, Chile; Pasquini et al. 2000). 
The fibers with UVES allowed us to obtain high-resolution spec- 
tra (R=47 000), for Red Giant Branch (RGB) stars, for red super- 
giant stars and for stars in the clump. The clusters were observed 
in Service mode in April, July, and August 2009, using a single 
configuration for each cluster The cross disperser CD3 was used 
resulting in a wavelength range ~4750-6800 A. A log of obser- 
vations (date, centre of field, date, observing time of each expo- 
sure, grating, number of exposures, number of stars) is given in 
Table |2] The spectra were reduced by ESO personnel using the 
dedicated pipeline. We analyzed the 1-d, wavelength-calibrated 
spectra using standard IRApQ packages and the SPECTRE pack- 



age by Sneden et al. (Fitzpatrick & Sneden [T987l ). One fiber for 
configuration was used to register the sky value, but the coiTec- 
tion was negligible. In Figure [1] we show a region of the spec- 
tra of stars in NGC 6192, NGC 6404, NGC 6583 (only mem- 
ber stars, see below). In Figures Fig. |2] Fig. [3] and Fig. |4] we 
present the color-magnitude diagrams of the three clusters. The 
stars observed spectroscopically are marked with filled (red) cir- 
cles. Stars which turn out to be non-members are marked also 
with a cross. 

First of all, we measured the Radial Velocities (RVs) using 
RVIDLINES on several tens of metallic lines on the individ- 
ual spectra. Using the RVs measured on the single lines, the 
task RVIDLINES produces an average radial velocity with its 
rms dispersion. We measured the radial velocity of each star on 
each individual spectrum, correcting them for the contribution 
of the different heliocentric velocities during each night of ob- 
servations. The uncertainty associated to the RVs of each star, 
corrected for the heliocentric component, is obtained from the 
rms value, after averaging estimates from different exposures; it 
turns out to be 0.5 km s"'. 

Finally we combined the spectra obtained during the differ- 
ent exposures with the task SCOMBINE. 

In Tables [3] ID |5j we show the available photometry of the 
observed stars, the signal to noise ratio (S/N), the RVs, corrected 
for the heliocentric velocity, and the membership derived on the 
basis of the RVs (M stands for "member", NM for "non mem- 
ber"). For NGC6192 the identification numbers and UBV mag- 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 



for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Table 2. Log of the observations. 



Cluster 


Centre of field 
RA Dec 
J2000.0 


Date 
yyyy-mm-dd 


Exptime 

(s) 


CD 


No. of exp. 
tot. 


No. of stars 


NGC6192 


16 40 50 


-43 22 00 


2009-04-09 


2775 


CD3 


1 


6 


NGC6404 


17 33 37 


-33 14 48 


2009-04-09 


2775 


CD3 


4 


5 








2009-04- i 1 


11 1 J 


CU3 












2009-08-06 


2775 


CD3 












2009-08-06 


2775 


CD3 






NGC6583 


18 15 48 


-22 08 00 


2009-04-09 


2775 


CD3 


5 


4 








2009-07-28 


2775 


CD3 












2009-08-20 


2775 


CD3 












2009-08-20 


2775 


CD3 












2009-08-20 


2775 


CD3 
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Fig. 2. Colour-magnitude diagram of NGC6192. The observed 
stars are marked with (red) filled circles. Stars which turn out to 
be non-members are marked also with a cross. The photometry 
is from Claria et al. (|2006] | 
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Fig. 3. Colour-magnitude diagram of NGC6404. Symbols as in 
Fig.|2] The photometry is from Carraro et al. (120051) . 



nitudes were taken from Clarfa et al. (2006'), for NGC6404 and 
NGC6583 the ID numbers and VI photometry are from Carraro 
et al. (I2005I I. For all clusters the JHK photometry is from the 
2MASS catalogue (Skrutskie et al. '2006V There are previous ra- 
dial velocity measurements for the evolved stars of NGC6192 
(Claria et al. 2006). We found good agreement with their val- 
ues. Their average velocity for five member stars, including two 
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Fig. 4. Colour-magnitude diagram of NGC6583. Symbols as in 
Fig.|2] The photometry is from Carraro et al. (I2005I I. 



binary stars (No. 91 and 96) is -7.7 + 0.38 km s"'. The aver- 
age velocity of our member sample is -8.9+0.7 km s"' in good 
agreement with the average velocity computed for the same stars 
using the Claria et al.'s measurements, -8.4+0.5 km s 

Previous radial velocity measurements are not available for 
NGC6404 and NGC6583. From our spectroscopy, the average 
velocities of member stars are 10.6+1.1 km s"' (4 stars) and - 
3.0+0.4 km s-i (2 stars) form NGC6404 and NGC6583, respec- 
tively. 

For our further analysis we thus exclude stars Nos. 4 and 
274 of NGC6192, No. 72 of NGC6404, and Nos. 12 and 82 of 
NGC6583 on the basis of their non-membership. 



4. Abundance analysis 

The method of analysis is very similar to that adopted in the pa- 
pers by Randich et al. (2006) and Sestito et al. (2006 20081), 
based on measured equivalent width (EWs) and on the use of 
the code MOOG. In this paper the analysis of chemical abun- 
dances was carried out with the version 2002 of the spectral 
program MOOG (Sneden 1973) and using model atmospheres 
by Kurucz (1993). Like all the commonly used spectral analy- 
sis codes, MOOG performs a local thermodynamic equilibrium 
(LTE) analysis. 
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NGC 6192 


Star 




DeC2MAS5- 


V 


B-V 


hlHASS 


HlMASS 


K2MASS 


RV 


S/N 


Notes 


















@6500A 


(km s-') 




0004 


16 40 45.63 


-43 24 07.8 


11.588 


1.334 


9.094 


8.345 


8.187 


170 


+46.3 


NM 


0009 


16 40 40.81 


-43 22 59.4 


11.430 


1.510 


8.478 


7.883 


7.645 


170 


-9.4 


M 


0045 


16 40 24.21 


-43 20 08.0 


11.743 


1.546 


8.674 


7.986 


7.797 


160 


-9.1 


M 


0096 


16 40 16.46 


-43 22 37.6 


11.306 


1.502 


8.277 


7.619 


7.371 


140 


-9.7 


M 


0137 


16 40 11.42 


-43 24 00.8 


11.223 


1.758 


7.776 


7.002 


6.757 


140 


-7.6 


M 


0274 


16 40 35.69 


-43 27 22.2 


11.300 


1.650 


7.948 


7.142 


6.888 


180 


+2.1 


NM 



Table 4. Data of observed stars in NGC6404. 



NGC 6404 


Star 


R-A2MA55 


DeC2M/lS5 


V 


V-I 






K2MA.V.V 


RV 


S/N 


Notes 


















@6500A 


(kms"') 




5 


17 39 30.86 


-33 14 52.3 


12.35 


2.379 


7.892 


7.139 


6.803 


120 


+9.7 


M 


16 


17 39 43.55 


-33 14 06.6 


13.42 


2.445 


8.716 


7.748 


7.427 


110 


+9.8 


M 


27 


17 39 43.46 


-33 15 14.0 


13.67 


2.433 


8.989 


8.061 


7.685 


150 


+ 10.3 


M 


40 


17 39 37.24 


-33 14 10.1 


14.09 


2.418 


9.471 


8.390 


8.098 


120 


+ 12.8 


M 


72 


17 39 34.92 


-33 14 39.3 


14.74 


2.293 


10.279 


9.385 


9.056 


100 


-42.9 


NM 



Table 5. Data of observed stars in NGC6583. 



NGC 6583 


Star 


RA2MA55 


DeC2M/lS5 


V 


V-I 


J2MA55 


H2MA55 


K2MASS 


RV 


S/N 


Notes 


















@6500A 


(kms-') 




12 


18 15 50.04 


-22 08 21.4 


13.05 


1.788 


9.663 


8.885 


8.663 


90 


+44.5 


NM 


46 


18 15 51.13 


-22 07 26.4 


14.10 


1.583 


10.986 


10.331 


10.170 


100 


-3.4 


M 


62 


18 15 51.23 


-22 08 28.0 


14.27 


1.652 


11.056 


10.379 


10.186 


80 


-2.6 


M 


82 


18 15 47.54 


-22 08 31.5 


14.49 


1.587 


11.105 


10.381 


10.198 


100 


+8.7 


NM 



4.1. Line list and equivalent widths 

Spectral lines to be used for the analysis of Fe I, Fe II and other 
elements (Mg, Al, Si, Ca, Ti, Cr, Ni) both in the Sun and in 
our giant stars were selected from the line list of Gratton et al. 
(2003) and Randich et al. (2006). The total list includes 161 and 
14 features for Fe I and Fe II, respectively, and several features 
of Al I, Mg I, Ca I, Si I, Ti I, Cr I, Ni I. In order to have a strictly 
differential analysis, we selected the same spectral range for the 
Sun and for our giant stars. 

We used the task eq in SPECTRE to normalize small portions 
of the spectrum (~40 A) and to interactively measure the EWs 
of the spectral lines by Gaussian fitting. Strong lines (EW > 150 
mA) have been discarded, since they could be saturated and they 
are critically sensitive to the microturbulence value, thus a more 
detailed treatment of damping would be necessary to fit their line 
wings. We also discard faint lines with EW < 20 mA because 
of the large uncertainties in their measurements. The values of 
EWs are available in electronic Tables A. 1-3, where the first two 
columns list the wavelengths and element, and the others show 
the corresponding EW for each star. 

4.2. Atomic parameters and solar analysis 

We adopted the oscillator strengths {log gf) of Sestito et al. 
(12006) . with exception of those of Al obtained by Randich et al. 
(I2OO6I 1 from inverse solar analysis. Radiative and Stark broad- 
ening are treated in a standard way in MOOG. We used the 



Unsold (11955 ^ approximation for collisional broadening since 
this choice does not affect the differential analysis with respect 
to the Sun as discussed by Paulson et al. ( I2003I I. In addition, very 
strong lines that are most affected by the treatment of damp- 
ing have been excluded from our analysis. The major differ- 
ence with the previous studies of Sestito and collaborators is in- 
deed in the different treatment of damping. Here we have chosen 
a more homogeneous treatment for all lines using the Unsold 
approximation with respect to the procedure by Sestito et al., 
who adopted two different types of damping coefficients, those 
of Barklem et al. (120001) . when available, otherwise the coeffi- 
cients of Unsold, multiplied by an enhancement factor In order 
to check the compatibility of the two approaches, we considered 
the EWs of two stars in NGC 6 253 (No. 023498 and 024707 
in Table 4 of Sestito et al. 120071 ) with stellar parameters similar 
to those of our stars. We recalculated [Fe/H] using our proce- 
dure, finding for No. 023498 and No. 024707 the [Fe/H] values 
+0.30 and +0.40, respectively. They compare well with the val- 
ues quoted in the original paper, i.e., +0.32 and +0.39, respec- 
tively. 

The first step is to derive the solar abundances of Fe and other 
elements, in order to fix the zero points of the abundance scale to 
minimize errors in the results. The solar spectrum was obtained 
with the same UVES setup of our observations. The line list for 
the Sun is available in electronic form (Table A. 4). The table 
includes wavelengths, name of the element, the excitation po- 
tential EP, log gf, and measured EW in the solar spectrum. We 
computed the elemental abundances for the Sun by adopting the 
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Table 6. Element abundances for the Sun in our analysis 
(Teff=5770, log^ = 4.44, ^=1.1), and in Anders & Grevesse 
(fTM9l AG89). 



El. 


No. 


logn(X)o 


log<X)o 




lin. 


Our 


AG89 


Mgi 


2 


7.56 ± 0.03 


7.58 


All 


2 


6.49 ± 0.03 


6.47 


Si I 


10 


7.58 ± 0.02 


7.55 


Cat 


9 


6.29 ± 0.03 


6.36 


Til 


14 


4.93 ± 0.03 


4.99 


Cri 


7 


5.67 ± 0.03 


5.67 


Fe I 


59 


7.47 ± 0.03 


7.52 


Fe n 


9 


7.49 ± 0.02 




Nil 


22 


6.24 ± 0.03 


6.28 



following effective temperature, surface gravity, and microturbu- 
lence velocity: Tefj - 5770 K, \ogg = 4.44, and ^ = 1.1 km s"'. 
In our analysis neither log n(Fe I) vs. EW nor log n(Fe I) vs. EP 
show any trend, implying that both the assumed microturbulence 
and effective temperature are correct, log n(Fe I) and log n(Fe II) 
are consistent within the errors indicating also a correct choice 
of the gravity value. Output solar abundances are listed in Table 
|6] together with those from Anders & Grevesse d 19891 1 that are 
used as input in MOOG. Our Fe I abundance is slightly lower 
than that measured by Anders & Grevesse (I1989I I. The origin 
of this discrepancy is probably due to the adopted log gf, but it 
would not affect the further analysis since it is strictly differen- 
tial. 

4.3. Stellar parameters 

We initially estimate effective temperatures (Teff) and gravities 
from photometry. We used the B-V and V-K vs. Teff calibra- 
tions by Alonso et al. ( fT999] l for the RGB stars of NGC6192, 
whereas we used V-K vs. T^ff calibrations by Levesque et 
al. (2006) for the post-RGB stars of NGC6404 and NGC6583. 
Surface gravities were derived using the expression 

log^ = log(M/Mo) + O.4(Mbol-Mbol0)+41og(reff/reff0) + log^O, 

where M is the mass and Mboi the bolometric magnitude 
(with MboiG = 4.72), Teff is derived as described above (with 
7^efFo=5770 K), and log ^©=4.44. The apparent magnitudes are 
transformed into absolute magnitudes using the distances and 
reddening values in Table 1. The bolometric corrections to be 
applied to My in order to obtain Mboi are taken from Alonso 
et al. (1999) in the case of giant stars and from Levesque et 
al. ( i2006) for post-RGB stars. The masses were derived from 
the isochrones of Girardi et al. ( 120001) using the estimated ages 
of each cluster (see Table |2l). They are: 4.1 Mg for NGC 6192 
(age ~ 0.18 Gyr), 3.0 Mq for NGC 6404 (age ~ 0.5 Gyr), 2.0 
Mo for NGC 6583 (age ~ 1 Gyr). 

The photometric Teff and logg values are the starting val- 
ues. For NGC6192, where two photometric estimates of T^s are 
available, we used as initial value their average. We optimized 
their values during the spectral analysis, employing the driver 
ABFIND in MOOG to compute Fe abundances for the stars: 
the final Teff estimate was chosen in order to eliminate possi- 
ble trends in [Fe/H] vs. ER The surface gravity was optimized 
by assuming the ionization equilibrium condition, i.e. log «(Fe 
II)=logn(Fe I). Then, if necessary T^s was re-adjusted in order 
to satisfy both the ionization and excitation equilibria. The other 
stellar parameter is the microturbulence ^ that was optimized by 




1 2 3 4 5 

E.P. 

Fig. 5. log «(Fe/H) vs. E.R for star No.40 in NGC 6404 with 
two Teff. 4060 K as derived from photometry and 4250 K as 
derived from the fulfillment of excitation equilibrium with the 
spectroscopic analysis with MOOG. The two continuous lines 
are the mean least-square fits to the data performed by MOOG. 



minimizing the slope of the relationship between Fe I/H and the 
observed EWs. 

The abundance values for each line were determined with 
the measured EWs and stellar parameters listed in Table IT] Final 
abundances for each star and each element were determined as 
the mean abundance from the different lines. Due to the large 
number of lines available, a 2cr clipping was performed for iron 
as the first step before the optimization of the stellar parame- 
ters. As akeady mentioned, [Fe/H] and [X/Fe] ratios for each 
star were determined differentially with respect to solar abun- 
dances listed in Table |6] 

The derived photometric and spectroscopic stellar parame- 
ters are shown in Tab. |7] col. 1 presents the star ID, cols. 2-4 
the photometric Teff and \ogg values, and cols. 5-7 the spectro- 
scopic Teft, log^ and ^ estimates. Finally, in cols. 8-9 the iron 
abundance and its dispersion are shown, whereas in col. 10 the 
number of Fe I and Fe II, in parenthesis, lines available for the 
spectroscopic analysis are presented. 

For NGC 6192 we found reasonable agreement between the 
photometric and spectroscopic parameters, while for NGC 6404 
and NGC 6583 the differences are important. The origin of this 
discrepancy derives from the large uncertainty on the reddening 
value and on the distance, due to the location of these clusters 
towards the GC. In fact, as far as T^ff estimates of these two 
clusters are concerned, we note that in order to satisfy the ex- 
citation equilibrium we had to consider Teff values larger than 
those derived from the photometric relations (see as an exam- 
ple Fig. |5j- These differences can be reduced by assuming that 
slightly higher values of the reddening with respect to those mea- 
sured by Carraro et al. (I2005I I through a photometric analysis. 
For both clusters, a value of AAy larger than the one quoted in 
Table 1 by 0.60 would reduce the difference among photomet- 
ric and spectroscopic Teff. The remaining differences in surface 
gravities can be attributed to several factors, such as errors in 
distance moduli. 
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Stai- 


Teff{B-V) 


Teff(V-K) 


log gphol 


f fspec 






[Fe/H] 


0-1 


No lines 




(K) 


(K) 




(K) 




km s ' 








NGC6192 


9 


5070 


5190 


2.34 


5050 


2.30 


1.75 


+0.19 


0.07 


59 (8) 


45 


5000 


5000 


2.38 


5020 


2.55 


1.60 


+0.08 


0.08 


119 (9) 


96 


5090 


5010 


2.23 


5050 


2.30 


2.10 


+0.13 


0.10 


77 (8) 


137 


4580 


4500 


1.90 


4670 


2.10 


1.80 


+0.07 


0.08 


91(9) 


Average [Fe/H] 














+0.12 


0.04 




NGC6404 


5 




4300 


1.96 


5000 


1.0 


2.60 


+0.05 


0.09 


63 (3) 


16 




4060 


2.08 


4450 


1.65 


2.10 


+0.07 


0.09 


67 (8) 


27 




4060 


2.17 


4400 


1.40 


1.80 


+0.20 


0.09 


47 (4) 


40 




4060 


2.29 


4250 


2.30 


1.40 


+0.11 


0.10 


89 (5) 


Average [Fe/H] 














+0.11 


0.04 




NGC6583 


46 




4625 


2.90 


5100 


2.95 


1.45 


+0.40 


0.12 


68 (8) 


62 




4510 


2.91 


5050 


2.75 


1.45 


+0.34 


0.12 


71 (10) 


Average [Fe/H] 














+0.37 


0.03 






!.5 6370 6370.5 
Wavelenght (A) 



6371 6371.5 6372 



Fig. 6. A region of the normalized spectra of member stars in 
NGC6404 centred in the [Fell] line at 6369.45 A. The continu- 
ous (red) spectrum is of star No. 5, which has a very low gravity, 
while the dotted line is the spectrum of star No. 45 in NGC 6192, 
having similar temperature and metallicity, but higher gravity. 



4.3.1 . Evolved stars in NGC 6404 and NGC 6583 

The criterion of selection of evolved stars to be spectroscopically 
investigated is based on their location in the color-magnitude di- 
agram. Due to different distance and reddening of the three clus- 
ters we did not select stars in the same evolutionary phase in the 
three OCs. In particular, the high reddening towards NGC6404 
and NGC6583 did not allow us to observe red clump giant 
stars, which lie at V>16. The evolved stars we have selected 
in NGC6404 and NGC6583 for our spectroscopic analysis (see 
Figs. [3] and |4]l are indeed in a latter evolutionary stage with 
respect to the RGB phase, moving towards the asymptotic gi- 
ant branch phase. The most extreme example is star No. 5 in 
NGC6404 (in the color-magnitude diagram of Fig. [3] this is the 
star with the lowest V magnitude). This star has an extremely 



low surface gravity combined with a high effective temperature. 
To show the reliability of our surface gravity measurement we 
plot in Fig.|6]a portion of the normalized spectra located around a 
Fe II line of star No. 5 of NGC 6404 and of a star with similar ef- 
fective temperature but with higher gravity. No. 45 in NGC 6 1 92. 
The EWs of Fe II lines are inversely related to the surface grav- 
ity, thus lower gravity implies larger EWs of Fe II. This is the 
case of star No. 5 showing very strong Fe II lines, for which we 
measured log g = I (see Table 7). Also the high microturbu- 
lence of star No.5 (2.6 km s ') is in agreement with its advanced 
evolutionary phase, being comparable with the microturbulence 
of red supergiant stars (~3 km s ', cf. Davies et al. i2009j . 



4.4. Error evaluation 



As usually done, we considered three sources of errors in our 
analysis: /) errors in measurements of the EWs; //) errors in 
atomic parameters; Hi) errors due to uncertainties in stellar pa- 
rameters. 

We estimated the errors in the following way. We consider 
the standard deviation (crl) around the mean abundance derived 
by individual lines as representative of the errors in EWs. We 
show the nns uncertainties in Table [T] for Fe I and Fe II and 
in Table |9] for the other elements. We consider as negligible the 
errors in atomic parameters because they are minimized by the 
differential analysis performed with respect to the Sun. We esti- 
mated internal errors due to uncertainties in stellar parameters 
(cr2) by varying each parameter separately, while leaving the 
other two unchanged. Following the error analysis by Randich 
et al. (2006), we assumed random uncertainties of +70 K, +0.15 
km s"', and ±0.25 dex in Teff, ^ and log g, respectively. We var- 
ied the stellar parameters of the warmest and coolest stars of our 
sample. No. 46 of NGC6583, and No. 40 of NGC6404, respec- 
tively. The variations in the abundance ratios due to changes in 
the stellar parameters are shown in Table [8] Then we assumed, 
as an estimate of the error due to the uncertainties on the stellar 
parameters, the maximum variation reported in Table[8]for each 
element. 
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Table 8. Random errors due to uncertainties in stellar parame- 
ters. 



No.46 


Tefl =5100 K, log g = 2.95,f = 


1.45 km s ' 


A 


ATeff = ±70 


Alogg = ±0.25 


A^ = ±0.15 




(K) 


dex 


1 -1 

km s 


[Fe I/H] 


0.04/-0.05 


0.02/-0.01 


-0.07/0.08 


[Mg/H] 


0.03/-0.03 


0.00/0.00 


-0.03/0.03 


[Al/H] 


0.05/-0.04 


0.00/0.00 


n m /n r\A 
-U.Uj/U.U4 


[Si/H] 


0.00/0.00 


0.05/-0.02 


r\ r\A if\ f\A 
-(J.(J4/U.U4 


[Ca/H] 


0.07/-0.07 


0.00/-0.05 


-(J.(Jo/(J.(J/ 


[Ti/H] 


0.09/-0.09 


0.00/0.00 


-0.07/0.08 


[Cr/H] 


0.08/-0.07 


0.00/0.00 


-0.07/0.08 


[Ni/H] 


0.03/-0.03 


0.04/-0.03 


-0.0 //0. 08 


No. 40: Tell = 4300 K, log g = 2.30, ^ = 


1.4 km s"' 


A 


ATeff = +70 


Alogg = ±0.25 


Af = ± 0.15 




(K) 


dex 




[Fe I/H] 


0.02/-0.02 


0.01/-0.01 


-0.05/0.07 


[Mg/H] 


0.01/-0.01 


0.05/-0.05 


-0.02/0.03 


[AI/H] 


0.02/-0.02 


0.01/-0.01 


-0.02/0.04 


[Si/H] 


0.07/-0.06 


0.07/-0.07 


-0.03/0.04 


[Ca/H] 


0.07/-0.06 


0.00/-0.01 


-0.08/0.08 


[Ti/H] 


0.09/-0.08 


0.01/-0.02 


-0.13/0.09 


[Cr/H] 


0.07/-0.05 


0.01/-0.01 


-0.08/0.10 


[Ni/H] 


0.03/-0.02 


0.07/-0.08 


-0.07/0.08 



5. Results 

5.1. Iron 

[Fe/H] abundances and rms values for the sample stars are listed 
in cols. 8 and 9 of Table |7] while in Fig. Q we show [Fe/H] 
as a function of the effective temperature. No notable trends of 
[Fe/H] as a function of Teff are present. The average metallicities 
for the three clusters are listed in Table Q They are all oversolar: 
[Fe/H]= +0.12 + 0.04 (NGC 6192), -i-O.ll ± 0.04 (NGC 6404), 
+0.37 ± 0.03 (NGC 6583). 

As our sample is made of OCs from the inner Galaxy, they 
are expected to be metal richer than the solar neighborhood, so 
that the estimated metallicities are in qualitative agreement with 
their Galactocentric distances. 



5.2. Other elements 

The [X/Fe] ratios for Mg, Al, Si, Ca, Ti (a-elements), Cr and Ni 
(iron-peak elements) are listed in Table |9] The number of lines 
available for each element in each star is shown in parenthesis. 
Due to the very small number of lines available for Mg and Al, 
we computed their ratios over iron line by line, and then we cal- 
culated the average. Errors on [X/Fe] values are the quadratic 
sum of the error crl for [Fe/H] and [X/H] values. 

Mean abundance ratios [X/Fe] together with 1-cr standard 
deviations are listed at the end of each group of stars in Table 
|9] No evident trends of [X/Fe] ratios with [Fe/H] are present. 
NGC 6404 shows a larger scatter than the other clusters due 
to the lower S/N of its spectra. The ratios of a-elements to Fe 
are close to solar in the sample clusters, with the exception 
of [Ca/Fe] in NGC 6192 that is slightly lower than the solar 

value, 0.2. Also the Fe-peak elements are solar, with [Cr/Fe] 

slightly below zero for NGC 6404 (—0.2). NGC 6404 shows a 
larger scatter than the other clusters due to the lower S/N of its 
spectra. 

In Fig. [8] we compare the a-elements Mg, Ca, Si, Ti over Fe 
ratios of our target clusters with the same ratios recently mea- 
sured by Bensby et al. (i2010i) for inner disk giants. The latter 



Q 5 1 NGC 65H3 
0,4 - 



0.3 - 



0,3 
1 0,2 



OJ 0,1 ^ 



1 _NGC 6192 



- 
4000 
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Fig. 7. [Fe/H] versus the effective temperature in the studied 
member stars of the three clusters: NGC 6583 (filled triangles), 
NGC 6404 (filled circles), NGC 6192 (filled squares). 



are located at Galactocentric distances of 4-7 kpc, hence con- 
sistent with the range spanned by our targets. The three clusters 
represent the metal rich extension of the mean trend of field gi- 
ants, however they do follow the trend quite smoothly. Only the 
most metal rich cluster, NGC 6583, might be perhaps a bit over- 
enhanced in Ca and Ti. We interpret this behaviour as an evi- 
dence that clusters and field stars belong to the same population, 
i.e., they have formed in the same, though extended, star for- 
mation episode. Note that the field stars in Bensby et al. (20ro]l 
are on average more evolved (mean log g = 1.5) than our clus- 
ter stars, therefore presumably slightly older, which is consistent 
with the higher metallicity of the (younger) cluster stars. 

6. The metallicity gradient of open clusters 

The radial metallicity gradient preserves the signature of the 
processes of galaxy formation and evolution. OCs are among 
the best Galactic stellar populations to study its shape and are 
unique tools to investigate its time evolution. For this purpose, 
a large sample of homogeneous data in terms of metallicity, 
age, and distance, is necessary. In this respect, our starting point 
is the database of OCs by M09, which is a collection of ele- 
mental abundances, ages and distances derived by several au- 
thors with high-resolution spectroscopy. In the database, the Rqc 
were taken from Friel et al. ( 1995 2002 2006 ) or calculated us- 
ing the cluster distance given in the original papers and a solar 
Galactocentric distance of 8.5 kpc, consistent with Friel et al. 
works. The ages for the clusters older than 0.5-0.6 Gyr were cal- 
culated in an homogeneous way using the morphological age 
indicator SV (Phelps et al. 1994) and the metallicity-dependent 
calibration of Salaris et al. (2004). For the younger clusters, M09 
used the most recent age determinations available in the litera- 
ture, such as those obtained with the lithium depletion boundary 
method. We updated the database of M09 with the latest high- 
resolution spectroscopic observations and with the results of the 
present work. The ages of the new clusters (only those older than 
~0.6 Gyr) have been recomputed homogeneously using the same 
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Table 9. [X/Fe] abundances and averages. Errors are the quadratic sum of crl on [X/H] and on [Fe/H]. 



Star 


[Mgl/Fe] 


[All/Fe] 


[Sil/Fe] 


[Cal/Fe] 


[Til/Fe] 


[Nil/Fe] 


[CrI/Fe] 


NGC6192 


9 

45 

96 

137 


-0.16+0.16(2) 
-0.01+0.16(2) 
-0.07+0.11(2) 
-0.04+0.11(2) 


-0.08+0.14(2) 
0.14+0.14(2) 
0.03+0.2(2) 
0.03+0.08(2) 


0.14+0.14(8) 
0.00+0.12(6) 
0.00+0.15(6) 
0.05+0.14(6) 


-0.08+0.13(11) 
0.02+0.12(11) 
-0.04+0.15(6) 
-0.09+0.13(5) 


-0.09+0.14(21) 
0.03+0.12(25) 
-0.02+0.13 (17) 
-0.04+0.12(16) 


0.02+0.13(34) 
-0.11+0.13(40) 
-0.05+0.13(34) 
-0.05+0.15(27) 


-0.09+0.14(10) 
0.04+0.11(12) 
-0.03+0.14(10) 
0.01+0.11(9) 


Average 


-0.07+0.05 


0.03+0.05 


0.05+0.05 


-0.05+0.04 


-0.03+0.03 


-0.05+0.03 


-0.02+0.04 


NGC6404 


5 

16 
27 
40 


0.10+0.10(1) 
-0.04+0.16(2) 
-0.05+0.16(2) 
0.00+0.11(2) 


0.05+0.20(2) 
-0.03+0.14(1) 
0.15+0.14(1) 
-0.05+0.08(2) 


0.23+0.15(5) 
-0.07+0.14(5) 
0.14+0.12(3) 
0.08+ 0.14(5) 


-0.24+0.15(3) 
-0.18+0.13(3) 
-0.26+0.12(4) 
-0.13+0.13(2) 


0.08+0.13(12) 
-0.15+0.14(11) 
-0.07+0.12(11) 
0.12+0.12(7) 


-0.05+0.13(16) 
0.13+0.13(10) 
0.17+0.13(14) 
0.00+0.15(29) 


0.20+0.14(4) 
-0.02+0.16(8) 
-0.02+0.11(8) 
-0.04+0.14(7) 


Average 


0.00+0.05 


0.03+0.07 


0.09+0.09 


-0.20+0.04 


-0.01+0.10 


0.06+0.08 


0.03+0.08 


NGC6583 


46 
62 


-0.08+0.16(2) 
-0.02+0.11(2) 


0.09+0.14(2) 
0.14+0.20(2) 


-0.05+0.14(3) 
0.06+0.15(3) 


0.00+0.13(7) 
-0.03+0.15(7) 


0.01+0.14(9) 
-0.03+0.13(12) 


0.02+0.13(28) 
0.10+0.13(22) 


-0.16+0.15(4) 
-0.15+0.14(7) 


Average 


-0.05+0.03 


0.11+0.02 


0.01+0.05 


-0.01+0.01 


-0.01+0.02 


0.06+0.04 


-0.15+0.01 



method adopted by M09. For most of the sample clusters (5V is 
available from the literature; when not available, we obtained 
this quantity from available color-magnitude diagrams. 

From the literature, we included 7 young clusters, namely 
NGC 6281 (0.3 Gyr), NGC 3532 (0.35 Gyr), IC2714 
(age~0.4 Gyr), NGC 2099 (0.4 Gyr), NGC 6633 (0.45 Gyr), 
NGC 1883 (0.65 Gyr), IC 4756 (0.79 Gyr), and 7 intermediate- 
age clusters NGC 5822 (1.16 Gyr), NGC 1817 (1.12 Gyr), 
Cr 110 (1.0 Gyr), Tombaugh 2 (2.13 Gyr), NGC 6939 
(2.05 Gyr), NGC 2158 (1.9 Gyr), and NGC 2420 (2.2 Gyr), 
and one old cluster Be 39 (7 Gyr). Abundance determina- 
tions are from Smiljanic et al. (|2U09] | (IC 2714, IC 4756 , 
NGC 3532, NGC 6281, NGC 6633), Jac obson et al. (l2009l l 
(NGC 1817, NGC 1883), Villanova et al. ( l20T0ll (T o 2), Pace 
et al. (gOlO) (NGC 5882), Pancino et al. OOTOl i (Cr 110, 
NGC 2099, NGC 2420), Friel et al. (20T0'(Be 39). W e also up- 
dated the M09 database with the results by Friel et al. (120101) for 
Be 31. NGC 6404 and NGC 6192 ai-e too young to be calibrated 
with the Salaris et al.'s relationship, while the age of NGC 6583, 
the oldest cluster of our sample, computed with the morpholog- 
ical age indicator 5V is ~1 Gyr, in agreement with the age given 
by Carraro et al. (2005) with isochrone fitting. 

6.1. The inner gradient 

The behaviour of the metallicity gradient in the inner part of the 
Galactic disk (Rgc <8 kpc) is still an open issue. The inner part 
of the disk is of paramount importance because it connects the 
chemistry of the disk with that of the bulge. 

Several studies identified one or more changes of slope in the 
[Fe/H] radial distribution. For OCs, Twarog et al. (1997) firstly 
identified a break in the [Fe/H] radial distribution of their cluster 
sample at Rgc ~10 kpc and this was recognized as the first in- 
dication of a transition between the inner and outer Milky Way 
disks (cf. Friel et al. 2010 ). The occurrence of a steepening of the 
metallicity gradient in the inner disk with respect to the outer re- 
gions was also suggested by Andrievsky et al. (120021 ) who found. 



using Cepheid abundances, two changes of slope in the radial 
gradient at Rgc ~ 11-12 kpc and Rgc ~ 7 kpc. The large 
compilation of Cepheid data by Pedicelli et al. (20091 ) shows in- 
stead a smooth and steady increase in the slope for Rgc < 8 kpc, 
without any further steepening towards the GC. They found an 
iron gradient in the inner disk (Rcc < 8 kpc) more than a fac- 
tor of three steeper (-0.130 + 0.015 dex kpc ') than in the outer 
disk (-0.042 + 0.004 dex kpc"'). Neither the metallicity gradi- 
ents of planetary nebulae (PNe), nor that of HII regions show 
any evident steepening in the inner regions (cf. Stanghellini & 
Havwood l20T0l Deharveng 200©. In the case of PNe and HII re- 
gions, gradients of oxygen and other a-element (i.e. S/H, Ne/H, 
Ar/H) were derived, rather than for iron; thus a direct compari- 
son is more difficult. In addition, few of these objects have been 
studied at Rgc > 15 kpc. 

In Fig.|9]we compare the iron gradient of OCs with that of 
Cepheids (data from Pedicelli et al. SOOp in the radial range 
4 kpc< Rgc <23 kpc. First of all, the metallicity scales of the 
two populations are in very good agreement. The gradients of 
young OCs and of Cepheids compare very well in the whole ra- 
dial range. A weighted mean least-square fit of the whole sample 
of young OCs (age<0.8 Gyr) gives for Rgc ^ 8 kpc a slope of 
-0.137+0.041 dex kpc-' andfor/^oc > 8 kpc-0.033+0.005 dex 
kpc"', which are in agreement with the values -0.130 + 0.015 
dexkpc ' and -0.042+0.004 dex kpc ' of Pedicelli et al. (20091). 
The main diff'erences between young (Cepheids and young OCs) 
and old (intermediate-age and old OCs) stellar populations re- 
side in the outer disk, where several old and intermediate-age 
clusters are observed. Their radial metallicity distribution be- 
comes flat at large Galactocentric radii, while this effect is not 
appreciable in the young populations due to the lack of mea- 
surements. 

Both populations show a clear steepening of the gradient at 
similar values of Rcc- However, due to the high dispersion at 
each Rcc> it is difficult to identify the exact radius where the 
change of slope happens. 



10 



Magrini, L. et al.: Chemistry and dynamics of OCs in the inner disk 




0.6 



0.2 



-0.2 



1 1 1 1 1 1 1 1 1 1 

• 

- • • 

•1 . •• • 

• M .. 

• • •• 


1 1 1 1 1 

• 
• 

^ * 


• 

: [Si/Fe] 

1 1 1 1 1 1 1 1 1 1 





-0.2 - 



0.6 




0.4 


• 






0.2 














— I 1 1 1 1 1 1 1 1 1 1 1 1 r 



V / • 

.*• 



[Ca/Fe] 



0.5 -1 
0.6 

0.4 

0.2 h 



-0.2 h 



-0.5 



0.5 



_l I I I I I I I I I I I I L_ 



T 1 — 1 1 1 1 1 1 1 1 1 1 r- 



••• 



• • • * 



[Ti/Fe] 



_l I I I I I I I I I I I I L_ 



-1 



-0.5 
[Fe/H] 



0.5 



1 



-0.5 
[Fe/H] 



0.5 



Fig. 8. The a- elements Mg, Ca, Si, Ti over Fe ratios of our target clusters (big blue filled pentagon) with the same ratios recently 
measured by Bensby et al. (120 lOt for inner disk giants (small red filled circles). 



As described in the Introduction, several chemical evolution 
models can reproduce changes of slope in the abundance gradi- 
ent (e.g., Lepine et al. 2003, M09, Colavitti et al. 2009 Fu et 
al.|2009) making different assumptions and hypothesis. Without 
entering in model details, the observational evidence of a steeper 
inner gradient suggests that the inner disk has evolved at a dif- 
ferent rate with respect to the outer disk. Many reasons can be 
at the origin of this; among others we recall the predominance 
of the halo collapse in the inner disk, the effect of the corotation 
resonance, the existence of density threshold in star formation, 
and/or radial varying star formation efficiency. Combinations of 
these causes are also likely to operate together. However, before 
a definitive interpretation in terms of models could be given, 
data still demand some further investigations and analysis. For 
example, there is still very little information about the possi- 
ble azimuthal variations of the gradient (e.g.. Luck et al. 120061 
Stanghellini & Havwood. l20T0l) . 

6.2. The time-evolution of the inner gradient 

In Fig. [To] we show the iron gradient in the very inner regions, 
between and 8 kpc form the GC. For Rqc < 4 kpc OC and 
Cepheid metallicities are not available, and we do not show 
metallicity measurements for the disk populations because of 



the difficulty to disentangle between disk and bulge stars. We 
compare the inner gradient traced by young populations, such 
as Cepheids and young clusters, with the gradient of older OCs. 
The [Fe/H] gradient of young clusters agrees with the Cepheid 
one, as expected for populations with similar ages and both 
NGC 6192 and NGC 6404 are in perfect agreement with respect 
to the average abundance of Cepheids. 

On the other hand, the [Fe/H] values of some old and 
intermediate-age OCs are surprisingly higher than for younger 
clusters; this is true in particular for the two innermost 
intermediate-age/old OCs: NGC 6583 (from this paper) and 
NGC 6253 (Sestito et al. l2007] l. The predictions of 'classical' 
evolution models for the Galactic disk, assuming that the disk it- 
self is formed by the collapse of the halo (with one or more pro- 
cesses of infall) and radial mixing is not considered, are in this 
respect well known. They imply that at each radius the metallic- 
ity increases with time, so that the older is the population con- 
sidered, the lower is the metallicity 

This is shown in Fig.[TO|in the case of the chemical evolu- 
tion model by M09 (see also, e.g., several cases of the two-infall 
model shown in Fig. 11 by Chiappini et al. 120011) . The M09's 
model reproduces only the disk evolution and does not include 
the bulge. The three curves represent the iron gradient predicted 
by the model at three different epochs in the Galaxy lifetime 
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Fig. 9. The iron gradient: OCs of different ages (top panel) and Cepheids (bottom panel). Empty circles (green) are OCs with 
age<0.8 Gyr, triangles (blue) with 0.8 Gyr<age<4 Gyr, and squares (magenta) with age>4 Gyr The three large filled circles (red) 
are the data presented in this paper. The black small empty circles in the lower panel are [Fe/H] of Cepheids by Pedicelli et al. 
(|2009] |. 



(present- time, 2 Gyr ago, 6 Gyr ago). For Rgc ^4 kpc the iron 
abundance decreases as the Galactocentric radius increases; at 
a given radius the metallicity at early epochs is lower than at 
present. Only in the very inner region, Rqc < 3 kpc, the gradient 
at early epochs and the gradient at present time are inverted, i.e. 
old population are slightly metal richer than young ones. This is 
due to the almost complete conversion of the gas into stars al- 
ready at early epochs, and to the subsequent dilution at recent 
times with gas expelled during the final phases of low-mass star 
evolution. 

One could imagine to change the model parameters in or- 
der to move the Galactocentric radius of the point of inversion 
of the metallicity gradient, and to reproduce the high metallic- 
ity of some of the old and intermediate-age clusters in the in- 
ner Galaxy. The model parameter that governs the position of 
the crossing-point is primarily the scale-length of the infall law 
(see M09 for detail). Varying the scale-length of the exponen- 
tial law driving the infall of gas from the halo, we can produce a 
more or less extended inner metallicity plateau. However reason- 
able values of the scale-length can vary the location only from 
~2 to ~4 kpc form the GC, and cannot reproduce the position. 



in the [Fe/H] versus Rgc plot, of clusters like NGC 6583 and 
NGC 6253. 

The location of these two clusters needs diff'erent explana- 
tions, such as, a different place of birth with respect to their 
present-time location or an incorrect (over-estimated) distance 
estimates. Non-circular orbits and/or radial migration of popula- 
tions (stars, but also gas) can explain the large scatter of [Fe/H] 
at each radius and at each age in the Galactic disk (see, e.g., 
the chemical evolution model of Schonrich & Binney 2009). In 
fact, not only clusters, but also field stars could have significantly 
non-circular orbits orbits, with apocentric and pericentric dis- 
tances differing by more than ~1 kpc. To verify this hypothesis, 
we have analyzed the orbits of the complete sample of inner disk 
OCs. 

7. The orbits of the inner OCs 

In this Section we describe the results of orbit calculations for 
a sample of inner-disk clusters. Firstly, we illustrate how abso- 
lute proper motion components are derived, and then we calcu- 
late clusters' orbits adopting a variety of Galactic potential mod- 
els. The main aim is to investigate whether, for clusters located 
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Fig. 10. The disk inner gradient (0 kpc < Rqc < 8 kpc). Symbols 
for OCs and Cepheids are as in Fig.|9] We show also the model 
curves of M09 at different epochs: present-time (continuous 
Une), 2 Gyr ago (dashed curve), 6 Gyr ago (dotted curve). 



inside the solar ring, the bar can significantly affect their orbit 
or not. For this result, among the various models, we will fo- 
cus only on the bar-less Allen & Santillan ( 1991 ) model and the 
FeiTers n - 2 model (Binney & Tremaine 2008), which has the 
maximum bar strength. Finally, we will use the orbital parame- 
ters to address the effect of the orbital motion on the chemical 
evolution of these clusters. Such an exercise is a classical one. 
Examples of that are in Carraro & Chiosi (1199 41). while open 
cluster orbits have been calculated several times in the past (see 
e.g. Bellini et al. [2010 or Carraro et al. 120061 for recent exam- 
ples). However, in all past cases, no effect of the bar could be 
seen, due to the mean large distance of the clusters from the 
GC and the relatively small sphere of influence of the bar. For 
the clusters under analysis here, no previous attempts have been 
done to derive their orbits. 



7.1. Proper motion components 

In order to obtain the clusters' mean absolute proper motion 
components and membership probabilities we took proper mo- 
tion data (and associated uncertainties) on individual stars from 
the UCAC3 catalogue (Zacharias et al. 1201 Oi l. Here we briefly 
recall the basics of the method. The interested reader can find all 
the details in Jilkova et al. (in preparation). 

The basic method to segregate cluster and field stars us- 
ing proper motion distributions was described long time ago by 
Sanders ( 119711 ). The fundamental assumption of this model is 
that proper motions are distributed according to circular or ellip- 
tical bi-variate normal distributions for cluster and for field stars, 
respectively. A maximum likelihood principle is then used to de- 
termine the distributions' parameters. Zhao & He dl 9901 1 later on 
improved the method by weighting proper motion components 
with observational uncertainties. 

We devised the 9 parametric models starting from the de- 
scription by Wu et al. (120021 ). but keeping the dispersion of the 
circular Gaussian distribution of clusters stars as zero. Therefore 



the width of the distribution is given only by the observed eiTors 
of proper motions (see also the discussion in Balaguer-Niinez et 
al.|2004). 

From the UCAC3 catalogue we selected only those stars 
having measurements of their proper motions and errors, and 
employed the 2MASS photometry as an additional check. The 
choice of the sample radius around each cluster center was made 
individually for each cluster taking into account the suggestions 
by Sanchez et al. (2010). To minimize the influence of possible 
outliers, we firstly considered all the stars belonging to one pop- 
ulation with a single bivariate normal distribution. Further we 
used only stars with both components of proper motion within 
the 3-cr interval. Then we restricted the range of the sample to 
\fi\ < 30 mas/yr. For such a sample the two-distributions analysis 
was done. The solution of maximum likelihood equations was 
found by a non-linear least squares minimization procedure. 

Once the membership probabilities from pure proper motion 
data were obtained, the CMDs using the 2MASS J magnitudes 
vs. (J-H) color indexes were inspected for the stars with clus- 
ter membership probability higher than 90%. For some clusters 
the additional magnitude constraint was used to eliminate field 
stars and iteratively the parameters of the new distribution were 
obtained. For more details and a comparison with previously de- 
rived UCAC2 proper motions, see Jilkova et al. (in preparation). 
In Table 10 we give the identification name of the clusters, their 
equatorial coordinates (RA, Dec), their Galactic coordinates (1 
and b), their distances from the Sun, their average radial ve- 
locities, proper motion in mas/yr (fj.a and fig), locations in the 
Galactic system (X, Y, Z) , radial velocity components (U, V, 
W), and finally age and metallicity [Fe/H]. 

7.2. Orbit computation 

With absolute proper motions, radial velocities and heliocentric 
distances of clusters we calculated their Galactic orbits. We as- 
sumed the heliocentric distances with 10% eiTors and we de- 
rived three sets of initial conditions for each cluster. The initial 
conditions were calculated following the procedure by Johnson 
& Soderblom (1987 ). We adopted the solar motion components 
with respect to LSR as bein g (U,V,W) = (11.1,12.4,7.25) km/s 
from Schonrich et al. (120101) . We also used a right-handed coor- 
dinate system where U is positive in the direction to the GC, the 
rotation velocity of the LSR of 220 km/s, and the Galactocentric 
distance of the Sun is 8.5 kpc. 

For the Galactic potential we adopted the axisymmetric time- 
independent model by Allen & Santillan (119911 ) and we also con- 
sidered an influence of the Galactic bar. For the latter we used 
a FeiTers potential (n = 2) of inhomogeneous triaxial ellipsoids 
(Pfenniger ll984l ). The characteristics of the bar model are: /) the 
bar length of 3.14 kpc; ii) axial ratios 10:3.75:2.56; Hi) a mass of 
9.0-10''^ Mq, iv) an angular speed of 60 km/s/kpc: and v) an ini- 
tial angle with respect to Sun-GC direction of 20°. We used the 
same parameters values as Pichardo et al. (2004) or Allen et al. 
(|2QQ8| for their more sophisticated models of the Galactic bar. 
For the axisymmetric background potential we kept the Allen 
& Santillan (1991) model, except for the mass of the bulge of 
4.26 X 10^ Modot (the bar replaces 70% of bulge mass). 

The integration routine uses a Runge-Kutta 4th order integra- 
tor (Press et al. 119921) . The relative change in the Jacobi constant 
is of the order of 10"^ to 10 We integrated the orbits back- 
wards in time over the intervals coiTesponding to the age of the 
clusters. For each orbit we give the most relevant orbital param- 
eters (mean peri- and apo-galacticon, maximum height from the 
plane, eccentricity, and birthplace) in Table [TT] Errors associ- 
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Fig. 11. Orbit computation for two clusters, NGC 3690 and 
NGC 6583, for the bar-less Allen & Santillan (1991) model 
(red), and the Ferrers n - 2 model (green). 



ated to each of these values are computed averaging the values 
obtained by running orbit calculations with 3 different distances. 
In FiglTT]we show, as an illustration, the orbit computation for 
two clusters, NGC 3690 and NGC 6583, for the bar-less Allen 
& Santillan (1991) model (red), and the Ferrers n - 2 model 
(green), to show a case where the bar does not have any influence 
(NGC 3690) and a case where the bar slightly affects the orbit 
(NGC 6583). 

In Fig [12] we show the efifect of the OC orbital motion on 
their radial position. For each cluster an horizontal bar shows 
the epicyclical variation from the apogalacticon to the perigalac- 
ticon. For some clusters, with eccentric orbits, large variations 
of their radial positions, up to about 4 kpc, are allowed during 
their lifetime. 

7.3. The influence of orbits on fhe shape of tlie gradient 

The most interesting quantity in order to derive the original 
shape of the metallicity gradient is the cluster birthplace. We 
have computed it integrating the orbits up to the epoch of birth 
of each cluster (see Table|8]l. The birthplaces, computed with the 
maximum age (Age + AAge) and with the minimum age (Age 
- AAge) are shown in Table [TT] We have assumed that AAge is 
~IQ%. The uncertainty due to the error on the cluster distance is 
estimated integrating the cluster's orbit with three different dis- 
tances, Rgc, Rgc ± A/?GC, and it is represented by the errors on 
birthplaces in Table [TT] 

In Fig. [T3] we show the metallicity gradient computed with 
three different sets of radial positions: /) the measured Rqc, H) 
the radial positions at the time of birth of the clusters from orbit 
calculations assuming the bar-less Allen & Santillan model, and 

the Rgc birthplace from orbit calculations with the n - 2 
Ferrers bar model. From Figs.[T2land[T3] we test our initial hy- 
pothesis of a different place of birth with respect to their present- 
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Fig. 12. The epicyclic variation of the radial position of the con- 
sidered OCs. In this plot we show the results from orbit calcula- 
tions for the bar-less Allen & Santillan model (upper panel) and 
for the bar Ferrers model (lower panel). 



time location of NGC 6253 and NGC 6583 and we check how 
much the present-time shape of the gradient is affected by cluster 
orbital motion. The initial radial position of some clusters dif- 
fers remarkably from their present-time position. In particular, 
for the clusters closer to the GC whose orbits might be affected 
by the presence of a bar, the choice of a model with or without 
bar is important to derive their birthplaces. 

NGC 6253 has a quite eccentric orbit and its radial position 
ranges from the perigalacticon ~5.4 kpc to the apogalacticon 
~8 kpc, both with bar and bar-less models (see also Fig. [TTT i. 
With an age estimate of ~4.5+0.5 Gyr, the birthplace of this clus- 
ter would be ~6 kpc with the bar model or ~7 kpc with the bar- 
less model, hence moved inward with respect to its present-time 
location. However, the age of NGC 6253 is particularly uncer- 
tain: Piatti et al. ( |1998i l found an age of 5+1.5 Gyr, Twarog et 
al. (^OOT) derived 3+0.5 Gyr, and Bragaglia & Tosi (2006) indi- 
cated an age of ~3 Gyr. The uncertainty on the age impacts on 
the position of the cluster along the orbit and might allow even 
more internal birthplace, up to ~5.4 kpc (the perigalacticon), as 
expected from its high metallicity. 

NGC 6583 is younger than NGC 6253, with an age of ~1 Gyr, 
thus the propagated uncertainty due to age determinations on 
its birthplace is smaller. In addition, the orbit of NGC 6583 is 
quite circular, with an epicyclic variation of ~1 kpc. With both 
Galactic potential models, the birthplace of NGC 6583 is ex- 
pected to have a small variation with respect to its present-time 
location. Thus NGC 6583 is really a high-metallicity cluster lo- 
cated at about Rgc =6.4 kpc. 

Weighted linear fits of the gradient slope are shown in 
Fig. [T3]for the Rcc. for birthplaces computed either with bar- 
less or with bar models. The gradient with measured Kcc has 
a steep slope driven by the position of NGC 6253. The three 
slopes are -0.177+0.034 dex kpc"', -0.110+0.026 dex kpc"', 
-0.1 14+0.021 dex kpc^', respectively. We excluded NGC 3690 
from the fits because its birthplace is in a different radial re- 
gion where the gradient changes its slope, becoming flatter. The 
inner gradient appears more defined and in better agreement 
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Fig. 13. The inner metallicity gradient of OCs with: measured Galactocentric distances (upper panel), radial position at the time of 
birth of the clusters, from orbit calculations for the bar-less Allen & Santillan model (middle panel), and from orbit calculations for 
the n=2 FeiTer bar (lower panel). The dotted lines are average least-squared fits. 



wi th Cep heids data (-0.130±0.015 dex kpc cf. Pedicelli et 
al. 120091 ) when we consider the initial radial position of each 
cluster, in particular for the birthplaces computed with the bar 
model. 

Thus, concluding, the knowledge of OC orbits is an impor- 
tant tool to derive correctly the shape and slope of the metal- 
licity gradient. Even if statistically the effect might be marginal 
(see also, e.g., Wu et al. 2009), the determination of clusters' or- 
bits and birthplaces is important to clarify the position of some 
'anomalous' clusters, such as NGC 6253, or to confirm their po- 
sition, as in the case of NGC 6583. 

8. Summary and conclusions 

In this paper we have presented new high-resolution spectro- 
scopic observations of evolved stars in three OCs located be- 
tween 6 and 7 kpc from the GC (NGC 6192, NGC 6404, 
NGC 6583) together with orbit computations for the complete 
sample of clusters with Rgc<8 kpc and available high-resolution 
spectroscopy. The UVES observations allow us to extend the 



sample of clusters studied in the inner Galactic disk and to dis- 
cuss the shape of the metallicity gradient in this radial region. 

The main results we have obtained ai^e the following: 

i) from radial velocity analysis, we obtained the membership 
of evolved stars finding 4 members in NGC 6192, 4 in 
NGC 6404, and 2 in NGC 6583; 

ii) from their high-resolution spectra and standard LTE analy- 
sis, we derived stellar parameters and abundance ratios of 
the iron-peak elements Fe, Ni, Cr, and of the a-elements 
Al, Mg, Si, Ti, Ca. Their average metallicities are [Fe/H] = 
+0.12+0.04(NGC6192), +0.1 1+0.04 (NGC 6404), +0.37 + 
0.03(NGC 6583) and [X/Fe] are consistent with solar values. 
In addition, we have compared our results with abundance 
ratios of inner disk giant stars finding a good agreement; 

iii) the clusters we have analysed, together with other OC and 
Cepheid data, confirm a steep gradient in the inner disk, a 
signature of a different evolutionary rate with respect to the 
outer disk. 
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iv) we found that, at a given radius, the metallicity of some 
intermediate-age and old clusters in this radial region is 
higher than that of young clusters (in particular this is true 
for NGC 6583 and NGC 6253). For the first time, we have 
coupled metallicity studies with dynamics, computing the or- 
bits of inner disk OCs and, from them, the birthplaces of 
the OCs with Rgc<8 kpc. We have found that there is a 
high probabihty that NGC 6253 had more internal value of 
Rgc than shown by its present position; this would recon- 
cile its high metallicity with its present-time radial location. 
NGC 6583 has instead a higher probability to be born close 
to its present-time position, and thus to be the most internal 
and metal-rich cluster cluster studied; 

v) the gradients computed using the cluster birthplaces instead 
of their present Rgc distances do not yield significantly dif- 
ferent values for the final slopes. However we found a bet- 
ter agreement with the Cepheid inner gradient slope and a 
higher correlation when we consider the birthplaces instead 
of the present-time Rgc- 
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NGC 6192 


16:40:23 


-43:22:00 


340.647 


2.122 


1500 


-7.70±0.38 


0.48±1.87 


1.43±1.67 


6.586 


-0.497 


0.056 


5.82 


237.15 


10.99 


0.20 


0.12±0.04 


NGC 6253 


16:59:05 


-52:42:30 


335.460 


-6.251 


1580 


-29.71±0.79 


-4.27±3.44 


-4.88±3.20 


6.571 


-0.652 


-0.172 


-36.75 


193.30 


13.12 


4.50 


0.36±0.07 


NGC 6281 


17:04:41 


-37:59:06 


347.731 


1.972 


480 


-5.58±0.26 


-2.30±2.11 


-4.57±1.84 


7.532 


-0.102 


0.016 


2.16 


215.26 


4.74 


0.30 
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NGC 6404 


17:39:37 
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-1. 177 
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-10.60±1.10 


0.75±2.80 


-2.32±3.48 


6.305 


-0.129 


-0.035 


-1.87 


213.53 


-1.19 


0.50 


0.11 ±0.05 


NGC 6583 


18:15:49 


-22:08:12 


9.283 


-2.534 


2100 


-3.00±0.40 


-1.12±2.68 


-1.45±2.76 


5.930 


0.338 


-0.093 


10.04 


206.99 


10.19 


1.00 


0.37±0.03 


NGC 6705 


18:51:05 


-06:16:12 


27.307 


-2.776 


1880 


35.08±0.32 


-4.98±2.86 


-3.05±4.13 


6.334 


0.860 


-0.091 


62.58 


202.57 


32.51 


0.20 


0.17±0.13 



Table 10. Parameters of the OCs located within ~8 kpc from the GC. 



Cluster 


Perigalacticon 


Apogalacticon 


Z,„™„ 


e 


Birthplace (age-Aage) 


Birthplace (age+Aage) 




kpc 


kpc 


kpc 




kpc 




Cr261 


4.87±0.23 


10.55+0.04 


0.161±0.011 


0.369+0.022 


7.20+0.63 


6.72±0.59 


1C4651 


7.21±0.08 


9.02±0.08 


0.168±0.004 


0.112+0.001 


8.45 ±0.48 


8.09 ±0.59 


NGC 3960 


6.62±0.14 


10.96±0.03 


0.248±0.013 


0.247±0.012 


9.61 ±0.30 


10.90 ±0.03 


NGC 5822 


7.02±0.11 


9.42±0.06 


0.050±0.002 


0.146±0.005 


8.62 ±0.50 


7.50 ±0.42 


NGC 6134 


6.94±0.11 


8.63±0.07 


0.069±0.008 


0.109±0.004 


8.15 ±0.38 


7.91 ±0.55 


NGC 6192 


7.07±0.08 


9.00±0.05 


0.128±0.003 


0.120±0.003 


7.64 ±0.00 


8.89 ±0.01 


NGC 6253 


5.37±0.19 


7.95±0.09 


0.210±0.007 


0.194±0.012 


7.51 ±0.23 


6.76 ±0.40 


NGC 6281 


8.03±0.03 


8.31 ±0.08 


0.055±0.001 


0.017+0.003 


8.25 ±0.08 


8.04 ±0.02 


NGC 6404 


6.76±0.11 


7.05+0.12 


0.078±0.007 


0.021+0.001 


7.01 ±0.11 


6.88 ±0.04 


NGC 6583 


5.91±0.15 


6.83±0.12 


0.131±0.004 


0.072±0.004 


6.28 ±0.21 


6.38 ±0.36 


NGC 6705 


4.91 ±0.1 8 


8.87±0.11 


0.364±0.018 


0.288±0.011 


8.74 ±0.03 


6.63 ±0.62 


Cr261 


4.86±0.23 


10.61±0.05 


0.261 ±0.0 11 


0.372+0.022 


6.36 ±0.07 


8.69 ±0.61 


IC4651 


7.21±0.08 


9.03±0.07 


0.168±0.004 


0.112±0.001 


8.49 ±0.46 


8.10 ±0.58 


NGC 3960 


6.64±0.13 


10.83±0.02 


0.247±0.012 


0.240±0.008 


9.25 ±0.53 


10.53 ±0.26 


NGC 5822 


7.02±0.11 


9.44±0.06 


0.050±0.002 


0.147±0.005 


8.66 ±0.52 


7.50 ±0.46 


NGC 6134 


6.92±0.11 


8.67±0.07 


0.069±0.008 


0.112±0.004 


8.22 ±0.34 


7.98 ±0.53 


NGC 6192 


7.09±0.08 


9.02±0.05 


0.128±0.003 


0.120±0.003 


7.62 ±0.01 


8.82 ±0.01 


NGC 6253 


5.34±0.19 


8.09±0.08 


0.211+0.007 


0.205+0.012 


5.43 ±0.16 


6.72 ±0.89 


NGC 6281 


8.03±0.03 


8.32±0.08 


0.055+0.001 


0.018+0.003 


8.26 ±0.08 


8.03 ±0.03 


NGC 6404 


6.68+0.13 


7.14±0.10 


0.078+0.007 


0.033+0.003 


7.09 ±0.05 


6.90 ±0.02 


NGC 6583 


5.86+0.21 


7.06+0.09 


0.132+0.004 


0.093+0.024 


6.21 ±0.28 


6.49 ±0.19 


NGC 6705 


4.92±0.17 


9.02±0.07 


0.366+0.018 


0.295+0.013 


8.75 ±0.06 


6.92 ±0.48 



Table 11. Output of the orbit calculation for the bar-less Allen & SantUlan model (upper section of the table) and the n=2 Ferrers bar model (lower section of the table). 
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